Thermodynamics second law and uj = —1 
crossing(s) in interacting holographic dark energy 

model 

H. Mohseni Sadjadi^*and M. Honardoost^ 
- - - ^ School of Physics, University College of Science, University of Tehran, 

■ North Karegar Ave. Tehran, Iran. 

. ^ Department of Physics, Shahid Bcheshti University, 

^ : Evin, Tehran 19839, Iran. 

X) 

^ : February 7, 2008 

O 

m 

J> ' Abstract 
^ . 

^ By the assumption that the thermodynamics second law is valid, we 

■ study the possibility of uj — —I crossing in an interacting holographic 

dark energy model. Depending on the choice of the horizon and the 
interaction, the transition from quintessence to phantom regime and 

I subsequently from phantom to quintessence phase may be possible. 

The second transition avoids the big rip singularity. We compute the 
dark energy density at transition time and show that by choosing ap- 
propriate parameters we can alleviate the coincidence problem. 
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1 Introduction 

X 

^ ' Recent observations suggest that the universe is undergoing an accelerated 

expansion [1]. This acceleration may be explained by the assumption that 
70% of the universe is filled by a perfect fluid with negative pressure, dubbed 
dark energy. Some present data seem to favor an evolving dark energy, cor- 
responding to an equation of state (EOS) parameter less than uj = —1 at 
present epoch from w > — 1 in the near past [2j. Many candidates for dark 
energy has been proposed such as the cosmological constant [3]: A constant 
quantum vacuum energy density which fills the space homogeneously, cor- 
responding to a fluid with a constant EOS parameter lo = —1; dynamical 
fields with a suitably chosen potential to make the vacuum energy vary with 
time [3], and so on. Recently, using holographic principle, a new candidate 
for dark energy which is independent of any specific field has been suggested 
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[5], [6]. Based on the holographic principle (which relates the number of de- 
grees of freedom of a physical system to the area of its boundary), in order 
to allow the formation of black holes in local quantum field theory, Cohen 
et al [7] proposed a relationship between UV and IR cutoff. This yields an 
upper bound on the zero-point energy density, which by a suitable choice of 
the infrared cutoff, can be viewed as the holographic dark energy density. In 
[3], three candidates for the infrared cutoff was proposed: the Hubble radius, 
the particle horizon and the future event horizon. There was shown that 
among these options only the future event horizon may be identified with the 
desired infrared cutoff. To study the coincidence problem and also to have 
other choices for the infrared cutoff, e.g. the Hubble radius, interaction 
between dark matter and dark energy may be considered in the holo- 
graphic dark energy model. As we have mentioned, based on astrophysical 
data, we may take into account the possibility of w = —1 (phantom divide 
line) crossing. Therefore dark energy models which can describe phantom 
divide line crossing, has been also studied vastly in the literature |9]. The 
phantom like behavior of interacting holographic dark energy was studied 
in [TU], where it was claimed that by selecting appropriate interaction pa- 
rameters the transition from the dark energy EOS parameter lod > —1 to 
ujd < — 1 is possible. Despite this, in [TT] it was shown that the dark energy 
effective EOS parameter cannot cross w^'^' = — 1. 

In this paper we consider interacting holographic dark energy model and 
study the ability of the model to describe the transition from quintessence 
to phantom regime and vice versa. After preliminaries in section two, where 
we introduce the interacting holographic dark energy model and some of 
its general properties used in the subsequent sections, in section three we 
study the possibility of crossing oj = —1. In section four we derive necessary 
conditions for existence of two transitions in our model. The first transition 
is from quintessence to phantom phase and the second the transition from 
phantom to quintessence regime. The importance of the second transition 
lies on the fact that it avoids the big rip singularity. We discuss also the 
behavior of Hubble parameter and dark energy density at transitions times. 

We use h = G = c = kB = l units throughout the paper. 

2 preliminaries 

We consider a spatially flat Friedmann-Lemaitre-Robertson- Walker (FRW) 
universe), with scale factor a{t) 

ds"^ = dt'^ -a^{t){dx'^ + dy^ + dz^). (1) 
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We assume that this universe is filled with dark energy and pressureless dark 
matter fluids satisfying the following equations of state 

P'd + 3iJ(l + ujd)pd = -Q 

p'm + 3Hpm = Q. (2) 

H = a/a is the Hubble parameter, pD is dark energy density and pm is 
the density of cold dark matter, "dot" denotes derivative with respect to 
the comoving time, w/j is dark energy EOS parameter. Q denotes the 
interaction of dark matter with dark energy. In this paper Q is assumed to 
be 

Q = {^mPm + \dPd)H, (3) 

where Am and \d are two real constants. Different choices such as Am = 0, 
Ad = and \d = Am has been adopted in literature [8], [10]. The Hubble 
parameter satisfies 

H = -Y\Pm + PD) 

and 

H = -Air ((1 + ujd)pd + Pm) 

= -47r(l+a;)p. (5) 

P = Pm + Pd is the total energy density satisfying 

p + 3H{l+Lo)p = 0, (6) 

and CO is the parameter of the EOS of the universe. In terms of ^la = Pd/p, 
we can write 

IVD^D = (7) 
Pd can be viewed holographic energy density 

The length scale L is an infrared cutoff and c > is a positive numerical 
constant. We assume 

L = PRfh + aRpH, (9) 
where Rfh and RpH are the future and particle event horizons 

f°° dt 

Rfh = a — 
Jt a 

dt 

RpH = a -, 10 
Jo « 
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and a and /? are two positive numerical constants. By taking /5 = and 
a = 1, we arrive at the holographic cosmology horizon adopted in [12]. For 
a = and /3 = 1, the infrared cutoff becomes the future event horizon [5]. 
Generally one can assume that L is a function of Kfh and Rpu [13]. The 
time derivative of L is obtained as 

L = HL + a- 13. (11) 

Using dH) and ([8]), HL in the above equation can be written as 

HL = c^}^K (12) 

By relating the entropy of the universe to the infrared cutoff L via 

S = ttL^, (13) 

the second law of thermodynamics results in L > 0, leading to pb < 0. In 
a phantom dominated universe (identified by i/ > 0), using ^ and ^ we 
get 

= ^ (pdH' - 2HHpd) < 0. (14) 
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For > we must have (LH) < which leads to L < 0. But if one 
requires L > 0, then either lim^^oo L = or L becomes positive after a finite 
time. 

In terms of the Hubble parameter, u is 

. = -1-^. (15) 

By substituting H = —H'^ + (/? — a)H/L + {HL)/L (which can be verified 
using (fTTl) ). ([TS]) becomes 

" = -3-3(^-")^ + 3m)^- ^''^ 



Using ([2]) one can show that 



r = 3rH 



I _|_ 1" 



(17) 



where we have defined r = Pm/PD- Therefore, by considering r = ^Ijj^ — 1, 
we obtain 

1 ^hD Ad Am 

'^^ ~ ~ 3H Qoil - ^d) ~ K'^ - ^d) ~ ^^d' 



and consequently 



1 Xpflp Am. . ^ 



4 



From (I16p and (I19p we can express oo and 17/) in terms of fi^): 



. = -1±^^^17, - ^(^4 - ^, (20) 
3 3c ^ 3 ' ^ ^ 



^ = (A™ - - l)f^?, + (1 - \^)^D + "^ 4(1 - ^i^)- (21) 

ii c 



For an accelerated universe, a{t) > (which results in a; < —1/3), we have 

(1 - A„) + (A^ -Xd- mo - ^ ^1) < 0. (22) 

In the absence of interaction, acceleration is not possible for /3 < a, as was 
pointed out by [5] for the case /3 = 0, a = 1. This was the motivation of 
[5] to take the future event horizon (instead of the particle horizon) as the 
infrared cutoff. However the above calculation reveals that in the presence 
of interaction, inflation may be possible even when the particle horizon is 
taken as the infrared cutoff. Combining (f2T]) and ([22]) . yields 

Qd < (23) 

Hence like the phantom regime, for /5 < a and in an accelerating universe, 
is decreasing. Applying the assumption L > 0, gives an upper bound 
for which depends only on H 

Qd < 2HVLd < 2H. (24) 

3 Crossing lj = —1 in interacting holographic dark 
energy model 

1 

At Lo = —1, u = $7 1) satisfies the cubic equation 

^^^^^u^ + {I + Xd - XmW + - 3 = 0- (25) 
c 

If w = — 1 is allowed, the above equation must have, at least, one positive 
root which is less than one. Based on Descartes rule, we know that the 
above equation, at most, has two real positive roots. So, if w = — 1 is 
crossed, two transitions may be possible, one from quintessence to phantom 
and the other from phantom to quintessence phase (by quintessence phase 
(regime) we mean u < —1/3 and H < 0). From p5|) it is clear that at 
a; = — 1 we have H = 0. If a transition from quintessence to phantom 
phase occurs at time ti, we must have H{ti) = 0, and H{t < ti) < 
and H[t > ti) > 0, therefore H{ti) must be a local minimum of H. In 
the same way, H must have a local maximum at t2, where the transition 
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from phantom to quintessence era occurs. In the neighborhood of ti, w is 
a decreasing function while in the neighborhood of ^ is an increasing 
function of time. So in order to see that \i uj = —1 crossing is permissible, 
we must also consider the behavior of the Hubble parameter near the roots 
of ([25]). 

In the following we assume that H > and L > 0. We also assume that 

the Hubble parameter is a differentiable function of time [H]. Therefore 
1 

following pT]) and p^ . L and u{= ^Ij^) are also differentiable. Let us 
consider the Taylor expansion oi H at t = ti, where ti is defined by H{ti) = 
(or uj{ti) = -1), 

H = hQ + hiT^ + 0{t''+^), a>2, (26) 

where r = t — ti, hi = ^^^i^{ti), and a is the order of the first nonzero 
derivative oi H a.t t = ti. If a is an even integer and hi > 0(< 0), then H 
has a minimum (maximum) at ti and the transition occurs at ti. Using (llSp 
we obtain 

o. = -l-||r'^-i + 0(0. (27) 
We consider the following expansion for u at t = ti, 



u = uo + uiT^ + U2-r''+^ + 0(r''+2), ni / 0, b>l (28) 

where b is the order of the first nonzero derivative of u at ti. If the solution 
(I26p is permissible, by inserting (I27p and (1280 in (I20p and by comparing the 
powers of r in both sides of (j20p . we obtain the following results: 

_ 1 + 2 _ ^(^^^3 = _i + ^. (29) 

^ ^ 3 ° 3c ° 3 ^ ' 

The roots of this equation specify u at transition time(s). In order that the 
transition occurs the above equation must have at least one real root in the 
interval (0, 1). (ii)For 

1 + Ap - Xm P-a , „ 

no / 0, (30) 



we obtain 



In the case 



ahi fl + XD-XmP-a^ , , ,„.^ 
uq uqUi, = a — 1. [61) 



3hl \ -3 



1 + Xd - Xm (3 - a 

5 Uq = 0, (32) 



and for (3 ^ a, we find 

2ahi {P — a)uQu1 



6h? 



= 26 + 1. (33) 

^ 
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In this case a is an odd number and transition does not occur. If /? = a, 
1 + Xd — Am = 0, and uj in (j20p . becomes a constant: uj{t) = —1, and no 
transition occurs. 

To determine a and the sign of hi, ([2T|) may be used : 



2u = H 



2{a-P) 4 , .X 3 2(/?-a) 2 X ^ 

^ + {\m - Ad - + ^-u^ + (1 - X^)i 



(34) 

By inserting (j28p and (|26p into the above equation, we see that the left hand 
side begins with t^~^, while if the right hand side does not begin with t^, it 
will begin by t^"/^''\ which is inconsistent with the left hand side. Thereby 
6=1. Hence based on ([^T|) and our previous discussion after (f^H|) . if the 
transition occurs, we must have a = 2. Note that for the case ()32p . we 
obtain a = 3 which may be corresponding to the inflection point of H. For 
6 = 1, by equalizing r'^'s coefficients in both sides of (pll) . we get 



2(a-/?) 4 , o 2(p-a) , ' 

^ ^ ' A + (Am - Ad - + ^ ng + (1 - Am)no 



2ui = /iQ 
which using ([29P reduces to 



(35) 



(3 — a 



no - 1 



(36) 



Now let us determine the sign of hi. Combining (I3ip and (I36p results in 



a/ii 
3^ 



1 + Ad - A, 



/? — Q 



-Mo 



^ — ^uo - 1 ) houl- (37) 



Prom (llip . (I12p and (I13p . by applying the thermodynamics second law we 
deduce 

l3 — a , , 

Mo - 1< 0, 38 

c 

which results in ui < 0. 5 = is ruled out because hi ^ 0. So we conclude 
that the sign of hi must be the same as the sign of (1 + Ad — Am)/(— 3) — 
(/3 — a)uo/c. For transition from quintessence to phantom phase we must 
have 

I + Ad-A^ P -a 

no > 0, (39) 

—3 c 

while a transition from phantom to quintessence era requires 

1 + Ad - Am P - a 



-no < 0. 



(40) 



Note that if, like |10j, we take a = 0, /? = 1 and Am = Ad, transition from 
quintessence to phantom is not allowed. 
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It is also instructive to study the behavior of L at tj. Assuming that 
-L > 0, by inserting the Taylor expansion of L at t = up to the order t"^ 
in (HI]), we obtain L = [Lq + Lit + O(r2)][/io + hir" + 0(r(''+i))] + a - /3, 
where Lq = L{ti) and Li = L{ti). Hence Li = h^L^ + a — (3. (fT2]) results in 
= c/{LH), therefore 

u — ^ u — ^ {\^^~^\ (41) 



which is consistent with (I36p . 

As a summary we have shown that in order that a transition phase 
occurs: (i) ()29p must have at least a positive real root in the interval (0,1) 
(ii) At these roots the Hubble parameter (if it is differentiable) must have the 
Taylor expansion ()26p . with an even integer a. In interacting holographic 
dark energy model we obtained a = 2 and verified that quintessence to 
phantom phase transition and vice versa occur provided ([5^ and ([in|) hold 
respectively. 

4 Two transitions in interacting holographic dark 
energy model 

In this part we try to investigate the ability of the system to return to the 
quintessence regime from the phantom phase. This may be interesting be- 
cause it avoids the big rip singularity which may be encountered in phantom 
models. 

us write ([29]) as 

Vul + gng + 1 = 0, (42) 



Let us write (|29[) as 
where 



2(/3-a) 1 + Ap - 
^ = T^^ = 5 — • 



In order to have two transitions, ()42p . must possess two real positive roots, 
which we denote by uqi, mo2i in the interval (0, 1). p and q are real numbers, 
hence the third root must be also real. Therefore the discriminant of 
i.e. — 27p^ — 4g^, must be positive 



p\ 2 / (7N 3 



2/ Mti <°- 

Prom X^j^j lioi^Oi = 0, we find that the third root, ■uqs, is negative and 
\uos\ < 1. So using < —U03U02U01 = 1/p < 1, we deduce p > 1. Also 
following Descartes rule of sign, having two positive roots is only possible 
when p > and q < 0. 
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The Sturm sequence constructed from (fl2]) is S{x) = {V{x),Vi{x),V2{x),V3{x)}, 
where 



V{x) = px^ + qx'^ + 1 
Vi{x) = 3px'^ + 2qx 

V2{x) = 



9p 

V^{x) = -^(27p2 + 4g3). (45) 



We have 5(0) = {1, 0, -1, -9p(4g3 _^ 27p2)/(4g4)|^ ^^^^ = {^ + ^ + 1^3^ + 
2g, -1 + 2g^/(9p), -9p(4g^ + 27p^)/(4g^)}. If we expect to have two real 
roots in the interval (0, 1), using ()44p and by applying the Sturm theorem, 
we find 

p + g + 1 > 

P Q 

- + - > 

2 3 



V3 



g\2 p 



3; i>»- 

These conditions and (I44p . may be resumed as 

3 

-<7-l<p<2(-|)% g<-3. (47) 

Now let us consider (|39p and (|4Up . A transition from quintessence to 
phantom phase may be occurred at ti {u{t = ti) = uoi), if /ii > 

l + Afl-Am /?-a 

UQi > 0. (48) 

—3 c 

For t > ti the system becomes phantom dominated until t = t2 {u{t = 12) = 
U02), i.e. when the second transition occurs, provided /ii < at t2 

l + Xn-^m 13 -a 

UQ2 < 0. (49) 

—3 c 

From (I14p we find that u is a decreasing function of time in phantom domi- 
nated era, hence uqi > U02 ■ So in order that (|18]) and (09]) become consistent 
we must have 

/3 - a < 0, 1 + XD-Xm>0- (50) 

The existence of uqi and uq2 which satisfy (|48p and (|49p . may be verified 
as follows: Using the condition posed on p and q, we obtain 
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The Sturm sequence at —2q/{3p) is 

S {-2q/{3p)) = {4<7V(27/) + 1, 0, -4q^/{27p^) - 1, -9p(V + 27/)/(4/)}. 

(52) 

So by invoking the Sturm theorem, one can verify that one of the positive 
roots, which based on our discussion in previous paragraph we take uo2, is 
located in (0, — 2g/(3p)) while the other, i.e. uqi belongs to (— 2g/(3p), 1). 
As we have previously mentioned (see discussion after (|25p ). w is a decreasing 
function of time in the neighborhood of uqi and an increasing function of 
time in the vicinity of uq2, whence for a differentiable uj, we must have 
a; = at a point in the phantom regime. It is interesting to note that 
for ti < t < t2, this happens at u = —2q/{3p), as it can be verified using 
([^U]) . Indeed from (pUj) . it is obvious that u;{t) = has three solutions: 
u = —2q/ (3p), u = and u = 0. In our case, u = is ruled out by ([23]) . The 
solution u = —2q/{3p), as we have seen lies in the phantom regime. The 
solution u = 0, must be obtained in the limit t — > oo in the quintessence 
era (by the assumption that the inflation will be continued). To verify this 
claim, we note that in an accelerating universe H + > 0, and for a > f3, 

(HL) = {H + H'^)L + {a- (3)H > 0. (53) 

Therefore HL is an analytic differentiable increasing function of time. Thereby 
u, is a decreasing function of time, as we showed before (see ([^H|) ) via another 
method. If n = occurs at a finite time t, then for t > t, u{t > t) < which 
conflicts with definition of u. In addition u[t) = leads to H{t)L{t) = oo 
which conflicts with continuity of H and L. 

As a test of our results, we have plotted a; + 1, ([20]) . in terms of u 
for an interacting holographic dark energy model with parameters {/3 = 0, 
a = 1, c = 1, Xd = 3.9, Am, = 2.5} (see fig.([l])). In this example w + 1 has 
two zero in the interval (0, 1), uqi = 0.86 corresponding to flo = 0.75 and 
uo2 = 0.72 corresponding to Qd = 0.52. For uq2 < u < uqi, the system is 
in phantom phase, i.e. lo < —1, and for u > uqi and u < no2 the system 
is in quintessence phase. Note that for an accelerating universe and when 
a > P, u is a decreasing function of time (see ([25]) and our discussion in the 
previous paragraph) and the directions of t and n-axis are opposite, uj = 
occurs in the phantom regime: u = 0.8, and at tt = 0. 

5 Summary 

In this paper we considered the holographic dark energy model with a general 
interaction between dark matter and dark energy, see We took the 
infrared cutoff as a linear combination of the future an particle horizon see 
® . We derived an expression for EOS parameter of the universe, uj, in terms 
of the ratio of dark energy density and total energy density of the flat FRW 
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Figure 1: cj + 1 as a function of u, for {/? = 0, a = 1, c = 1, Ad = 3.9, 
Am = 2.5} 

space-time, Qd, see (f20|) . Using ([20]) and the expression obtained for time 
derivative of fi/j (I2ip . and by assumption that the thermodynamics second 
law is stiU vahd, we studied the possibihty (and necessary conditions) for 
quintessence to phantom phase transition and vice versa with a differentiable 
Hubble parameter. Using some theorems about solutions of cubic equation 
satisfied by 0^ see (|25]) . we showed that such transitions occur provided 
we appropriately choose the parameters of the system see (j^7|) and ([^0]) . 
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